Dilepton flow and deconfinement phase transition in heavy ion collisions 
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The dilepton radial flow in Au+Au collisions at ^/snn —200 GeV is investigated. The space-time 
evolution of the fireball is described by a 2 + 1 dimensional ideal hydrodynamics with a variety of 
equations of state. The slope parameters of the transverse momentum spectra from the partonic 
and hadronic phases show distinct features and are sensitive to equation of state parameters. The 
elliptic flow and breaking of Mt scaling are also studied and have distinct features for the two 
phases. These features can serve as clean signals for the formation of a quark-gluon plasma in 
ultra-relativistic heavy ion collisions. 



I. INTRODUCTION 



Among all observables for determining the quark gluon plasma (QGP) in heavy ion collisions, the electromagnetic 
probes such as photons and dileptons are expected to provide clean signitures due to their instant emissions once 
produced [H-tllj. These thermal photons and dileptons contain undistorted information about the space-time trace 
i— i of the new state of matter formed in such collisions. There are many sources of dileptons in heavy ion collisions. 
' In the lower invariant mass region (M <1 GeV) dileptons are mainly from resonance decays and may be related to 
I \ chiral symmetry restoration |12l4l5j . In the higher invariant mass region {M > 3 GeV) dileptons are dominated by 
q ■ the Drell-Yan process and charmonium decays. For moderate invariant mass dileptons (1 < M < 3 GeV), it was 
^ \ argued that dileptons from semileptonic decays of correlated open charm in pp collisions are dominant [l6| . But in 
Ph . Au+Au collisions, both charm related single lepton contributions and their dynamic correlations are expected to be 
suppressed by medium modification. Therefore thermal radiation may play important role in the intermediate mass 
. region and the dilepton spectra can be used to extract thermodynamic parameters of the fireball. 
J> ' The observation of jet quenching and strong elliptic flow at the Rclativistic Heavy Ion Collider (RHIC) at 
\ Brookhaven National Laboratory tell us that the dense matter produced at RHIC interacts strongly and may reach 
local thermalization in very short time (l7l [l8j , implying that ideal hydrodynamic models are applicable to such 
systems [l9|, [2(| ■ In contrast to hadronic flow, the observables of dileptons are more direct and penetrating probes to 
. i the early space-time profile of the QGP [t| [l(| • The radial flow of thermal dileptons has been measured at CERN 
super-proton synchrotron (SPS) by the NA60 collaboration [2l|. It was found that the inverse slope parameter T e // 
increases with the invariant mass M of the lepton pair below the p meson mass and then starts decreasing above 
M ~ 1 [2l[ . The reason for the drop of T e f f around M ~ 1 GeV is not fully understood although it was thought to be 



. . , an indication of the transition to an emission source with much smaller flow possibly a partonic source [2ll . |22j . The 
^ ■ observed strong correlation of T e f f versus M in the region M < 1 GeV is mainly due to the collectivity developed 
k>( \ in the hadronic stage at SPS energies. Should the collectivity have been developed in the partonic phase of the 
evolution in high energy nuclear collisions, RHIC and/or LHC, one would also expect to see the increase of T e ff in 
the intermediate mass region 1 < M < 3 GeV. Inspired by the NA60 result and in order to develope clean observables 
for the formation of the QGP, we propose to study the transverse momentum distributions of di-clectrons over the 
entire region of 0.5 < M < 3 GeV. 

In this paper we use a 2+1 dimensional ideal hydrodynamic model to give the space-time evolution of the medium 
created in Au+Au collisions at ^/syvw =200 GeV. Our program gives results consistent with AZHYDRO (23|. To 
include the pre-equilibrium emission of the di-electrons at very beginning, we set the initial time for the hydrodynamic 
evolution tq = 0.2 fm/c as in Ref. Q instead of t$ = 0.6 fm/c in previous studies. The initial transverse energy 
density is calculated in the Glauber model with the peak temperature being at about 520 MeV in central collisions. 
Another critical input is the equations of state (EOS) of the dense matter (24|. We choose four types of EOS, (i) 
QGP-EOS, the ideal gas EOS for 3-flavor QGP, e = 3n = (19/12)tt 2 T 4 without the hadronic phase; (ii) HG-EOS, the 
resonance hadron gas EOS for the hadron gas (HG) [25[ without the partonic phase; (hi) MIX-EOS, the one with the 
first order phase transition with both the partonic and hadronic phases [25j ; (iv) LAT-EOS, the one extracted from 
the lattice calculations j26|. Note that in the MIX-EOS and LAT-EOS there are partonic and hadronic components. 
These EOS are shown in Fig. [T] We put emphasis on the lattice EOS which can describe the crossover from resonance 
HG to QGP in the temperature range 180-200 MeV. Other EOS are also considered in our calculations for the purpose 
of comparison. Althrough the crossover does not have a rigorous critical temperature T c to separate the QGP from 
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FIG. 1: Equations of state: QGP-EOS (upper-left), HG-EOS (upper-right), MIX-EOS (lower-left), LAT-EOS (lower-right). 
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the HG phase, we still choose T c = 180 MeV as a tuning parameter. The freezeout temperature is set to Tf = 120 
MeV below which there is no di-electron emission from the HG sector. The post-freezeout di-electron decay of the 
p mesons is not included directly in our calculation, whose effect on the di-electron spectra can be partly taken into 
account by lower the freezeout temperature. The fine tuning of T c and Tf does not qualitatively change our results 
and conclusions. 



II. DILEPTON YIELDS 

In the thcrmalizcd medium HG or QGP, the rate for the di-lepton production per unit volume is given by 



diN a I ( 2m?\ /. 4m? _ _ fl , . 

n B(P ■ «) 1 + -rA V 1 - T75- Imn (?)• (!) 



d 4 xd 4 p 4tt 4 M 2 ^ ' \ M 2 J V M 

Here m; is the lepton mass, a = e 2 /(Ait) is the fine structure constant with the electric charge e for leptons, p = 
(po,p) =Pi+P2 is the dilcpton 4-momcntum and M = \[p 2 , tib{p-u) — \ j {e v ' u l T — 1) is the Bose distribution function 
with T and u the local temperature and medium velocity respectively, is the retarded photon polarization tensor 
from the quark- loop or the hadronic loop, and H R = |n^. For the hadron gas, is further related to the retarded 
rho-meson propagator D R via Indl^ = —{e 2 rn A p /g 2 p )l l a\D R 1 where g p is the photon-rho-meson coupling constant in 
the vector meson dominance model, and m p is the rho-meson mass. The imaginary part of the retarded p-meson 
propagator is given by, 



ImlP? 

r m n fi = p (2\ 

p (p 2 -m 2 p + Rcn«) 2 + (Imn«) 2 ' 1 j 



where is the retarded rho-meson polarization tensor contraction. We assume that the hadron gas is mainly 
composed of mesons and consider following mesonic vertices for the p meson @ in n«: pirX and pK #'(1270), where 
X denotes meson resonances below 1300 MeV, namely, u, /ii(1170), ai(1260) and 7r'(1300). Here we neglect the 
baryonic contributions since the collisional energy considered in this paper is the RHIC energy (200 GeV) where the 
baryon yields are smaller than meson yields. The baryonic contributions may broaden the width of the rho-meson 
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even more. But we argue that they will not qualitatively change the behavior of T e f f in the mass range considered in 
this paper, considering that the mesonic resonance contributions do not qualitatively change the dilepton spectra and 
then the behavior of T e ff compared to the dominant irir process in this mass range. Note that 11^ depends on thermal 
distributions of mesons and then on local temperature T(x) and medium velocity u(x) determined by hydrodynamic 
simulation, which are functions of space-time position. We only take into account the processes mediated through the p 
meson, because the di-electron production is dominated by the isovector channel instead of the isoscalar one. According 
to the flavor SU(3) quark model, the relative weight of the electromagnetic couplings for vector mesons V = PjU>, 4> is 
about 9:1:2, roughly consistent to the electromagnetic decay widths Ty^ ee = 7.0,0.6, 1.27 KcV respectively On 
the other hand, w and <fi mesons can be easily identified in experiments. We also neglect the Daliz decay channels for 
77 and 7T°: 77 — > e + e~7, /i + /i~7 and ir° — > e + e~j. The pion spectrum is mainly below m, = 135 MeV and irrelevant 
to our current range of M. The 77 contribution can be easily deducted as the background in experiments due to its 
very long lifetime (about 1.5 x 10 5 fm/c), leading to decays outside the freezeout scope. 

After making approximation ns(p • u) w e p ' u ^ T and then integrating over the lepton pair rapidity y, we obtain the 
differential cross section for an expanding fireball (27| . 



d*N a 1 2m 2 \ I Am 2 



p T dp T MdMd(j) p 2tt 4 M 2 V M 2 J V M 2 
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where px = \Pt\ is the scalar transverse momentum of the lepton pair, <j) p and <f> v are the azimuthal angles of 
and the local transverse fluid velocity Vj- respectively, 7t = 1/ yl-t^ is the local transverse Lorentz factor for 
the fluid element, and Kq is the modified Bessel function of the second kind. Note that the encoding of space-time 
history of the fireball is realized by integrals over fluid coordinates d A x = rdrdr/d 2 ^, where r, 77 are proper time and 
space-time rapidity respectively and is the transverse position of the fluid element. For central collisions with 
azimuthal symmetry, the angular integral can be worked out by J d<fi exp [a; cos <f>] = 2ttIq(x) analytically with Jo being 
the modified Bessel function of the first kind. 

The evolution of the rate with the proper time r is shown in the left panel of Fig. [2] The emission rate is 
proportional to T 4 St where S is the transverse area. For the harder EOS like the QGP ideal gas there is a rapid 
expansion leading to a very fast decreasing of T and S. As a result the QGP freezes out in a very short time with less 
dilepton emission. But with softer EOS like the lattice one, T decreases slowly as S expands, making the freezeout 
isothermal lines expand and remain almost constant at large radii, which greatly increases the di-electron emission in 
the HG phase. 

The numerical results for the invariant mass spectra of the multiplicity with the LAT-EOS in central collisions are 
shown in the right panel of Fig. [2] The contributions from the QGP and the HG components arc distinguished. 
The shapes of the spectra with the QGP-EOS only or the HG-EOS only are similar to the QGP or HG component 
(bluc-dashed/red-dottcd line) here. One can see that the QGP/HG contribution dominates in the region M ^ 1 GcV. 
The emission from the QGP phase in the early stage with high temperatures contributes to the hard parts of the 
spectra with large M and tot- The HG phase in later stage with low temperature contributes in small M region and 
is shaped by the p meson form factor, but a sizable collective flow developed in this stage can harden the transverse 
momentum spectra [27i ] . Also we have shown the HG contribution from the prnr vertex only, which is lower than the 
full HG contribution (with all resonances below 1300 MeV included) in the lower mass region. This indicates the 
broadening effect to the p meson spectra from the resonances' contribution. 



III. EFFECTIVE TEMPERATURE, ELLIPTIC FLOW AND M T SCALING 

To discuss the sensitivity of dilepton signals to different stages and collectivity developed in the history of the 
fireball, it is valuable to investigate dilepton transverse momentum spectra and its inverse slope parameter. We can 
use the reduced transverse mass mj- = Mr — M with My = \J M 2 + p\, to replace pr as the variable. After carrying 
out the space-time integrals we assume that the transverse spectra can be approximately parameterized as [HI, [29| , 



d 2 N lTV^TM, ew (_^ + M l (4) 
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FIG. 2: (Color online) Differential multiplicity as a function of proper time(left panel) and the di-electron invariant mass (right 
panel). The LAT-EOS [2(| is used in the calculation and the contributions from the QGP and HG are shown in dashed and 
dotted lines, respectively. The dashed-dotted-dotted line is the HG contribution from the prnr vertex only. The contribution 
from the HG is dominated by p mesons. 




FIG. 3: (Color online) Left panel: slope parameter T e // as a function of M with the LAT-EOS and for three mr values. Results 
for three tut values are shown as solid, dashed and dotted lines. For mr = 0.5 GeV, thin red and blue lines with circles and 
triangles represent the results extracted from contributions of the HG and QGP phases, respectively. Right panel: The slope 
parameter T e // calculated for four equations of state at fixed uit = 2.5 GeV. 




with the average temperature T, the average transverse velocity vt of the fluid and its Lorentz factor 7 T = 1/ \J\ — v\. 
The asymptotic forms of the slope parameter T e ff can be written as, 

j T + M*v T , for p T <M 

where M* is a monotonic function of M. As emissions in different space-time are encoded in the spectra (j4|), the 
slope parameter T e ff depends on the fitting window of M and pt, so do T and vt- 

The slope parameters T e ff versus M from the LAT-EOS are shown in the left panel of Fig. [3j For low px, e.g. 
iriT = 0.5 GeV, T e ff approximately follows T e // ~ T + M*v\ as in Eq. ([5]). The curves below 1 GeV have smaller T 
and larger vt, while above 1 GeV they have an opposite trend. For even higher M the increase of T e // is not caused 
by the collective flow since larger M reflects the higher temperature in earlier emission when the collective flow has not 
fully developed. The disconnected curves with solid circiles and triangles are T e ff extracted from the contributions 
of the HG or QGP phase only. They reflect different trends of T e ff below/above 1 GeV from the HG/QGP phase. 
When rriT is larger, e.g. = 2.5 GeV, T e tt does not follow the above formula for low p?, but the same trends 

still exist. For very large mj, e.g. tot = 5.5 GeV, T e // follows the blue-shift formula T e // ~ T^J jjjfjjj -, which is 
independent of M. The valley at M ~ m p is due to the fact that the QGP component (with higher T e tf) is dominant 



5 



FIG. 4: Slope parameter T e // as a function of M with the dependence on the initial time to (left panel) and on the phase 
transition temperature T c (right panel). We choose rxiT = 0.5 GeV and use LAT-EOS for both panels. The red and blue lines 
with circles represent the results extracted from contributions of the HG and QGP phases, respectively. 




over the HG one (with lower T e //) except in the region near M ~ m p . One can see in the figure that T e ff increases 
with rriT, since the larger tut probe the earlier state of the fireball with high temperatures. For the MIX-EOS, since 
there is a large contribution from the coexisting stage of two phases with the same temperature and fluid velocity, the 
difference in T e fj between the QGP and HG phases is smaller than the case of the LAT-EOS. For other two types 
of EOS, the QGP and HG ones, T e // simply increase with M monotonously. In right panel of Fig. [3j we compare 
T e ff with different EOS. It is found that the magnitude of T e ff for the HG-EOS is much smaller than that for other 
EOS with partonic phase. Such behaviors of di-electron T e ff are expected to be measured and tested in the future 
experiments. 

In Fig. |4]we show the dependence of T e ff on the initial time to for hydrodynamic evolution and the phase transition 
temperature T c . When we tune To from 0.2 fm/c to 0.6 fm/c and change the initial peak temperature from 520 MeV 
to 350 MeV, or in other words, delay the hydrodynamic process, the magnitudes of T e ff, T and vt decrease due to 
smaller initial temperature and acceleration but their structure remains the same. When changing T c from 180 to 
150 MeV, the magnitude of T e // does not change much, but the slope increases in the small M region. This can be 
understood because a smaller transition temperature means that the HG phase occupies outer layer of the fireball 
with larger radial flow velocities. 

Now we can compare our result for T e f / with NA60 data [2l[ . The T e j f curve of NA60 is extracted from low nir 
spectra. It is comparable to our result in Fig. |4]with tot = 0.5 GeV. But the trend above M = 1 GeV in NA60 data 
is different: it is downward while our curve is upward after a transition region between 0.8 and 1.2 GeV. We now give 
the reason for such a difference. We note that the region of M > 1 GeV probes earlier stage of the fireball, where the 
QGP emission dominates over the hadronic emission. Such a different trend is just the manifestation of the radial 
flow and average temperature effects of the QGP phase. At the NA60 energy, the life of the QGP phase (if there is 
any) is much shorter and its average temperature is lower (near freeze-out Tf), the radial flow does not have enough 
time to develope and the average temperature is just above Tf . For M < 1 GeV, the dominant contribution is from 
the hadronic phase which appears at later stage of the evolution and has larger collective flow, so T e f f is larger since 
it has a component proportional to the flow velocity squared times an effective mass, see first line of Eq. ([5]). That 
is the reason for the upward/downward trend below/above 1 GeV at the NA60 energy. In contrast, at the RHIC 
energy, the life of the QGP phase is longer and there is some time for the radial flow to develope and make T e ft 
larger. Another effect is that the average temperature in the QGP phase is larger than at the NA60 energy. This the 
reason for the upward trend above 1 GeV in our result. 

For non-central or peripheral collisions we can compute the elliptic flow coefficient vi versus transverse momentum 
and invariant mass, V2{pt,M) = (cos(2</> p )). Setting the impact parameter 6 = 7 fm, we calculat dilepton V2 as a 
function of invariant mass for four EOS, see the left panel of Fig. [SJ The results with the MIX-EOS or LAT-EOS 
have similar features to Ref. 0. Dilepton sources within different mass windows arc sensitive to different stages 
of the expansion history. We see that t>2 is much smaller for larger M than that for smaller M. This is because 
dilcptons at larger M come mainly from the QGP phase at earlier time of the fluid evolution with smaller collective 
flow, while those at smaller M are mainly from the HG phase at later stage when the collective flow is fully developed. 
The HG contribution is characterized by the peak at p meson mass which quantitatively different from the partonic 
contribution. If there is only one phase in the course of the evolution, these distinct feature will disappear. The 
initial temperature of the QGP phase in EOS containing the QGP phase is much higher than that of the pure HG 
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FIG. 5: Elliptic flow «2 (left panel) and the ratio R for the Mt scaling (right panel) for different EOS. For «a the impact 
parameter is set to b — 7 fm. For R we scan all values of the impact parameter, i.e. b — [0, 15] fm. 
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phase in HG-EOS at the same initial energy density. So the early stage emission from the former is much larger than 
that from the latter (since yield is proportional to T 4 ), then the later stage emission for HG-EOS has much larger 
proportion in the total yield. This helps us understand the bigger elliptic flow of HG-EOS than that of other EOS 
containing the QGP phase in the intermediate mass region. As the temperature window for emission with HG-EOS is 
not very wide and the low mass dileptons come mostly from the low temperature stage, similar to the surface emission 
of hadrons, the elliptic flow has a slight increase with mass in the low mass region. The peaks at the p meson mass 
in the MIX-EOS and LAT-EOS cases and the increase in the small mass region are partially caused by the surface 
emission, but the main reason is the influence of the QGP phase in the region below the p meson mass [j| . 

We also note that if the QGP is created, the dilepton spectrum in the mass region between the </> meson mass and 
J/ip mass depends approximately only on its transverse mass Mt- This can be seen from Eq. ([3]), after the integration 
over azimuthal angle <j> p , the differential yield is proportional to Iq^tVtPt /T)Kq{^tMt /T). In the continuum region, 
the lepton mass, the quark mass and the collectivity are negligible, then the yield is only a function of Mt- This 
property is called the Mt scaling [3Cj| . It can be broken when transverse expansion becomes appreciable or hardonic 
gas emission with extra mass scales dominate the yield. By looking at the ratio, 

R _ dN/dM T dp T dy\M T =2.6GeV, PT =2GeV ^ 
dN/dM T dp T dy\ M T=2.6GeV,p T =0GeV ' 

we study the Mt scaling with more realistic hydrodynamic evolution. As shown in the right panel of Fig. [5l the 
scaling is rather robust if the QGP phase is present. We find that the elliptic ex pan sion, the collective flow of the 
QGP and the p meson form factor do not qualitatively change the result in Ref. [3(J], if the di-electron yields from 
partonic/hadronic source dominate in the intermediate/low mass region respectively. 



IV. SUMMARY AND CONCLUSION 



In conclusion, we have investigated the thermal production of di-electrons in the quark-gluon plasma created in 
heavy ion collisions at RHIC energy. Different types of the equations of state are used in our study. The phase 
transition from hadron gas to quark-gluon plasma leads to a rich structure for the thermal dilepton production. The 
mass dependence of the inverse slope parameter is sensitive to the collective dynamics of the medium: at the mass 
region (M <1 GeV), it is dominated by the hadronic interactions while in the intermediate mass region (1 < M < 3 
GeV), partonic interactions become important. At lower transverse momenta, the slope parameter for mixed-phase 
and lattice EOS shows different trends in M in the low mass region (hadronic phase) and intermediate mass regions 
(partonic phase), which reflects the existence of two distinct phases. In this case, there is transition area around 
M = 1 GeV to connect two T e ff components with different trends. In the hadronic phase in the lower mass region 
the flow velocity is found to be stronger than that in the partonic phase in the higher mass region. Around M ~ 1 
GeV, the slope parameter is found sensitive to the equation of state of the fireball used in our calculation. The elliptic 
flow and the Mt scaling are calculated which show distinct features for the hadronic and partonic phases. We have 
investigated the surface emission or the emission in a fixed interval of temperatures, where we have found that these 
features are absent for the two phases. In this sense the dilepton emission is a volume effect which can carry the 
information of the space-time history of the fireball. Therefore all these features of the collective flow can serve as 
clean probes to hot and dense medium created in high energy heavy ion collisions. 
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